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ABSTRACT 

We have studied the X-ray emission within the inner ~150 pc radius of NGC 
4151 by constructing high spatial resolution emission line images of blended OVII, 
OVIII, and NelX. These maps show extended structures that are spatially corre- 
lated with the radio outflow and optical [OIII] emission. We find strong evidence 
for jet-gas cloud interaction, including morphological correspondences with re- 
gions of X-ray enhancement, peaks of near-infrared [Fell] emission, and optical 
clouds. In these regions, moreover, we find evidence of elevated NelX/OVII ra- 
tios; the X-ray emission of these regions also exceeds that expected from nuclear 
photoionization. Spectral fitting reveals the presence of a collisionally ionized 
component. The thermal energy of the hot gas suggests that > 0.1% of the 
estimated jet power is deposited into the host interstellar medium through inter- 
action between the radio jet and the dense medium of the circum-nuclear region. 
We find possible pressure equilibrium between the collisionally ionized hot gas 
and the photoionized line-emitting cool clouds. We also obtain constraints on 
the extended iron and silicon fluorescent emission. Both lines are spatially un- 
resolved. The upper limit on the contribution of an extended emission region to 
the Fe Ka emission is < 5% of the total, in disagreement with a previous claim 
that 65% of the Fe Ka emission originates in the extended narrow line region. 
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Introduction 



Multi- wavelength high spatial resolution observations of NGC 4151 (e.g..lHutchings et al. 
I999I : lYang et al.lboOll : iMundell et aPboOsI : IWang et al.lboOQbi Istorchi-Bergmann et aliboool ) 
have provided unique information on the contin uum and emission lin e morphology in the 
central few 100 pc {d ~ 13.3 Mpc, 1" = 65 pc, iMundell et al.lll999| ). contributing to our 
und erstand i ng of the nature of the nuclear emission and its interaction with the host galaxy 



see 



UlrichI l200d for a review). In this paper we investigate further these topics, using 



high quality deep Chandra ACTS observations of the nearby prototypical Seyfert 1 galaxy 
NGC 4151, which allow the morphological investigation of spectral line emission in the inner 
circum-nuclear region (r ~ 150 pc in projection). In particular, we will address the radio 
jet-host galaxy interaction and the nature of the Fe Ka line emission. 

High resolution radio surveys show broad existence of collimated rad io outflows or jets in 



Seyfert galaxies, albeit less prominent than those of radio galaxies (e .g. , lUlvestad fc Wilson 
19841 : iKukula et al.lll995l : lNagar et al.lll999l : IXanthopoulos et al.boiol ). Ever since the Hubble 
Space Telescope emission line imaging r esolved the NLR mo r phology in a numb e r of Seyfert 
galax ies with linear radio structures (jCapetti et al.l 1 19951 : iFalcke et al.l Il998l : ICecil et al. 



20001 ). it has been actively debated whether the radio jet plays a competin g role against 



nuclear photoionization i n the narrow-line region (NLR) ionization structure ( Whittle et al. 



1988 



1999 



Rossi et al 



Binette et al.lll996l:lDopita fc Suther 



200 



3 



Whittle fc Wilson 



2004 



and 1996:lBicknell et al.lll998l:IWilson fc Ravmond 



Rosario 



20071 ). C/ian(ira X-ray observations 



of the nuclear regions of nearby Seyfert galaxies appear to offer powerful diagnostics for 
discriminati ng between photoionization by the n uclear radiation or coUisional ionization by 
a radio jet (ISako et al.ll2000l : lYoung et al.ll200ll ). althoug h the importance of radio jets in 



shaping NLR environments seems to differ ca se -bv-case (iWhittle et al.l 120051 : lEvans et al. 
20061 : ICandhi et allbood : iRosario et allboiOal lbl: iBianchi et allboioh . 



The nuclear region of NGC 4151 i s known to host a two - sided, ^300 pc-long radio outflow 



Tr 

along position angle (P.A.) ~77° (e.g.. lJohnston et al.lll982l:lPedlar et al.lll993t IMundell et al. 



I995I ). Within the central 100 pc of the galaxy, IMundell et al.l (120031 ) identifled a faint, 
highly collimated jet underlying discrete components that are shock-like features associated 
with the jet-gas clouds interactions. Although the distribution of the ionized gas and the 
kinematics in the NLR imply that shock ionization originated from jet-cloud i nteractions is 



unlikely to be the dominant source of ionizing photons compared to the nucleus (iKaiser et al. 
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2000l : iMundell et all 120031 ). Mundell et al.l (120031 ) find evidence that some NLR clouds are 
responsible for bending the je t. Near-IR emission-line mapping of the nuclear region by 
Storchi-Bergmann et al.l (|2009[ ) also reveals enhancement of [Fell] emission, suggesting the 



prese nce of shock heating in addition to nuclear photoionization ( IStorchi-Bergmann et al. 
2009h . 



In our previous work using the Chandra HRC image of NGC 4151. IWang et al.l (j2009b( ) 
found an excess of X-ray emission at these jet-cloud interaction locations. However, no 
spectral information could be obtained because of the very limited energy resolution of 



(Ogle et al. 


2000; 


Yang et al. 


2001) 



did not allow detailed high spatial resolution comparison to the radio jet and emission line 
gas, nor spectral studies of X-ray emission from spatially resolved f eatures. Our 200 ks ACIS- 



S observation of NGC 4151 (PI Fabbiano; see IWang et al.ll2010al . and Wang et al. 2011 for 



previous work on this data set) allows us to pursue these investigations, by comparing the 
morphology of the NelX and OVII with that of the radio jet. 

Moreover, with these data we can address the outstanding controversy on the spatial 
properties of the 6.4 keV Fe Ka emission of NGC 4151. In Seyfert galaxies, this emission 
line is expected to originate from cold m atter near the nucleus (< 1 pc), either the obscur- 



George fc Fabian 



1991 



ing torus (e.g., iKrolik fc Kallmarull987f ) or an accretion disk (e.g. 

Revnoldslll997l ). These model predictions, however, conflicted with the finding of lOgle et al. 
(I2OOOI ). who reported spatially resolved narrow iron Ka line emission in the Chandra HETG 
observation, concluding that 65%±9% of this emission originates in the ENLR at distance up 
to 6" from the nucleus (~400 pc acr oss). This con c lusion has been more recently contested 
by the ZMM- Newton based work of ISchurch et all (120031 ). 



2. Observations and Data Analysis 



The details of our ACIS observations and data reduction are described in IWang et al. 



(I2OIII . hereafter Paper I). Briefly, NGC 4151 was observed by Chandra for a total of 180 ks 
(after screening for high background in tervals) with the spec troscopic array of the Advanced 
CCD Imaging Spectrometer (ACIS-S; iGarmire et al.l l2003l ) in 1/8 sub-array mode during 
March 27-29, 2008. The data were reprocessed following standard procedures, using CIAO 
(Version 4.2) with the CALDB 4.2.1 provided by the ChandraX-iSij Center (CXC). Subpixel 
event repositioning ( "static" method in Li et al. 2004, using the corner split and 2-pixel split 



^See the Chandra ABC Guide to Pileup available at 



http : / / CXC . harvard . edu/ ciao/ download/ doc/pileup_abc . pdf 
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events) and subpixel binning techniques were applied to the ACIS images to improve the 
spatial resolution. 



During our previous work ( iWang et al.ll2010bl ) and Paper I, we became well aware of 



the complexity in data analysis caused by the bright nuclear emission of NGC 4151. First, 
photon pile-up is present in the nuclear region even with the reduced frame-time of our 
observation. We established that for the soft X-ray emission, pile-up is moderate (< 10%) at 
r > 1". Second, although the Chandra!s PSF is highly centrally peaked, the contamination 
from the scattered nuclear emission to the extended emission is not negligible due to the 
brightness of the nucleus and the broader PSF in the higher energy range. We performed 
Chandra PSF simulations that provide an estimate of the expected contamination from the 
nuclear emission in an extended feature. Detailed explanation of the analysis leading to 
these conclusions is presented in Paper I. We have taken into account this information in 
the following analysis. 



2.1. Soft X-ray Emission Line Images 



The 2-7 keV emission is dominated by the unresolved nucleus ( Yang et al.ll200ll : IWang et al. 



2009b|). Since we are interested in the extended X-ray emission, we used the energy range be- 
low 2 keV where this emission is prominent (see Paper I). To investigate the general spectral 
dependency of the morphology of this soft X-ray emission, images in three spectral bands 
below 2 keV were extracted from the merged data: 0.3-0.7 keV ("soft b a nd"), .7-1.0 keV 
("medium band"), and 1-2 keV ("hard band"). Following IWang et al.l (l2009a[ l. exposure 



maps were created for the individual bands to obtain exposure-corrected flux images. 

Figured] presents the resulting false color composite image of the central 7" x 7" (450 pc 
on a side) region of NGC 4151, where the soft, medium, and hard band images are shown in 
red, green, and blue, respectively. The images in the three band have been smoothed with 
a FWHM = 0.3" Gaussian kernel. It clearly shows bright structured soft X-ray emission 
along the northeast (NE) - southwest (SW) direction, which is also the direction of the radio 
jet. In particular, the medium band image shows a jet-like ~2" linear extension (contours 
in Figured]). This medium band is dominated by the NelX emission (see below). In the 
following we concentrate on the analysis of this inner 2"-radius circum-nuclear region. 



High spectral resolution grating spectra of NGC 4151 (lOgle et al.ll2000t ISchurch et al. 



2004J : lArmentrout et al.l 120071 ) have shown that the emission in these soft spectral bands is 
almost entirely dominated by lines. Although the spectral resolution of ACIS CCD cannot 
provide unique identifications of the strongest soft X-ray emission lines (<2 keV) seen in 
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by the HETG observations (I Ogle et al 



the nuclear spectrum (IWang et al.ll2010bl). we can identify the dominant transitions guided 



20001 ) ■ Most notably the blended lines appear as 



three strong lines in the ACIS spectrum, approximately c entered at . 57 keV , 0.68 keV, 
0.91 keV (c.f. Figure 4 in Yang et al. 2001; Figure 3 in IWang et al.1 boiObl ). We then 
extracted the X-ray emission in three narrow energy intervals (0.53-0.63 keV, blended OVII 
f,i,r; 0.63-0.73 keV, blended OVIII Lja and OVII RRC; 0.85-0.95 keV, blended NelX f,i,r) 
to create line strength images, highlighting regions of these prominent emission lines. This is 
a reasonable approach since the line em ission dominates over the weak underlying continuum 
( lOgle et al.ll2000l : ISchurch et al.ll2004 ) in these i iarrow bands and reveal substructures that 



are n ot obvious in the broadband images (e.g., 
2010h . 



Bianchi et al. 2010: Gonzalez-Martin et al. 



The resulting images are shown in Figures |2] and [31 showing the X-ray emission-line 
structure close to the nucleus (the central r = 2", ~130 pc) in the context of the radio 
outflow and the optical emission line clouds. The alignment between these images is done 
using the peak of X-ray emission, [OIII] emission, and the radio core. Since the astrometry 
of our Chandra image is accurate to 0.3", our interpretation of these features is not affected 
by significant alignment uncertainties. The position of the nucleus is indicated with a cross. 
The OVII, OVIII, and NelX emission line images all show extended morphology and some 
structures, clo sely following the P.A. of the large sc ale extended NLR traced by the [OIII] 



emission (e.g. 



Winge et al.lll997l : iKaiser et al.ll2000[ ) 



In particular, the linear feature seen in the medium-band is clearly present in the 
NelX image. There appears to be two X-ray enhancements bracketed by the optical clouds 
(indicated by arrows in Figure [3]) that are close to the radio knot features C2 and C5 
(iMundell et al.lll995l ). These "hot spots" are bette r visualized in the N elX/OVII ratio im- 
age, shown in Figure H together with the VL BA jet (IMundell et al.ll2003[ ) and the NIR [Fell] 
line emission ( IStorchi-Bergmann et al.ll2009l ). The locations where the jet appears to be in- 
tercepted by the optical clouds and where prominent [Fell] emission arise, are characterised 
by a NelX/OVII ratio of 2.8 ± 0.2, which is significantly (~ Qa) higher than the ratio in 
the surrounding regions (1.2 ± 0.2). We also examined the optical clouds with high velocity 
dispersion, possibly associated with jet-cloud impact (marked as crosses in Figure IH see 
Figure 7 in Mundell et al. 2003), and find an elevated NelX/OVII ratio of 2.1 ± 0.4 (with 
marginal significance). The eastern NelX/OVII hot spot appears to be associated with clear 
interactions between the jet and the gas cloud, where shock heating becomes important. We 
note that the observed line flux of OVII is more suppressed than NelX in the pr esence of the 
Galactic absorption column {Nh = 2 x 10^° cm~^), and XSPEC (jArnaudlll996l ) simulations 
show the absorption corrected ratio may systematically decrease by 10%. Thus the presence 
of NelX/OVII hot spots is not affected, unless there is significant differential absorption in 
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the region. 



2.2. Spectral Analysis of the Soft Emission 



We extracted X-ray spectra from both the high NelX/OVII emission line ratio regions 
and the surrounding low ratio region (Figure HJd), and fit them jointly (with the same model 
components but normalizations set free) for comparison using XSPEC (Version 12.5; Arnaud 
1996). Spectra and instrument responses were generated using CIAO tool spece^raS The 
background spectrum is taken from a nearby source-free region on the same CCD node. 
Spectra were grouped to have a minimum of 20 counts per energy bin to allow for fitting. 
The contribution from the bright nuclear emission cannot be neglected here; to have a self- 
consistent model, we included a fraction of the scattered emission predicted from the PSF 
simulation of the point-like NGC 4151 nucleus. 



by 

i 



We have made use of the Cloudy photoionization modeling code, last described b 
Ferland et al.l ( 119981 ) . to model the soft X- ray emission . Usi ng Cloudy version C08.0 
which enables a Cloudy/XSPEC interface (jPorter et al.l l2006l ) . we attempted to produce 
the soft part of the X-ray spectrum assu ming an open plane-paralle l geometry ("slab"). 
The dimensionless ionization parameter ( lOsterbrock fc Ferlandl |2006| ) is defined a.s U = 
Q /{AT[r'^cnH)i where uh is the hydrogen number density, r is the distance to the inner face 
of the model slab, c is the speed of light, and Q = j^^^y Ly/hu is the emitting rate of 
hydrogen ionizing photons (s~^) by the ionizing source. It was previously noted that the 
spectral fit to the soft X-ray continu um of the nucleus is quite uncertain because of heavy 
absorption (lArmentrout et al.l 120071 ) , therefore we adopted the broken power- law form in 
Kraemer et al.l ( l2005l ) for the AGN continuum {Gamma = 2.3 for the energy range between 
13.6 eV and 0.5 keV, and Gam ma = 1.5 for E > 0.5 keV). Normalizing to the observed 
X-ray luminosity of the nucleus (IWang et al.ll2010bl ). we obtain Q ~ 7 x 10^^ (photons s~^). 
We varied U and the column density Nh of the model slab to create spectral model grids, 
which were fed to XSPEC. 

We quickly find that, using a single photoionized component alone cannot reproduce 
the observed soft X-ray emission in NGC 4151. With a high photoionization component 
(log U ~ 1.9), we were able to produce the hydrogen-like neon and oxygen lines em ission (e.g.. 



OVIII, NeX), which are unambiguously present in all the grating spectra (e.g., lOgle et al. 



^http : //cxc . harvard . edu/ciao/threa ds/specextract7| 
^Available at http : //www . nublado . org/wiki/DownLoad| 
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2000l : ISchurch et al.l 120041 : lArmentrout et al.l 120071 ). However, the helium-like transitions, in 



particular the OVII and NelX emission line features, are not fitted accepta bly. Adding a sec- 
ond lo wer ionization photoionized component (log?7 ~ 0), as detected in lArmentrout et al. 
(120071 ). we were able to obtain a significantly improved fit (Ax = 200), except for residuals 
in the 0.7-1.1 keV range that are particularly strong for the hot spot spectrum. Adding 
additional photoionized components to the model does not further improve the fit when an 
F-test is performed. Allowing different absorption columns for the models also does not 
give a better fit (reduced xl, = 1-7), indicating that the high NelX/OVII ratio is not due to 
higher obscuration. The absorption column require d by the fits is consistent wit h the low 
Galactic extinction towards the NLR region found in ICrenshaw fc Kraemerl (l2005n. Iiistead , 
these residuals disappear when a thermal emission component [APEC; ISmith et al.ll200l[ ) 
with a temperature of kT = 0.58 ± 0.05 keV is added to the model, suggesting the presence 
of coUisionally ionized Fe L emission. The spectral fitting results are summarized in Table 1 
and shown in Figure [5l The absor ption column is consist e nt wi th the low Galactic extinction 
towards the NLR region found in I Crenshaw fc Kraemerl ( 20051 ). 

Lastly, we note that fitting the data with only combinations of absorbed optically-thin 
thermal emission (APEC) gave strong residuals of emission line features (reduced ^ 1), 
even when the abundance Z is allowed to vary. 



2.3. [OIII]/Soft X-ray Ratios in the Substructures 



The subpixel resolution ACTS image also allows a similar comparison of the X-ray emis- 
sion with the optical NLR clouds resolved in the HST image. Taking advantage of the 
spectral cap ability of ACIS, we measured the 0.3-2 keV counts for the same [OIII] clouds 
identified in lWang et al.l ( l2010bl ). and derived their 0.5-2 keV X-ray fiuxes using the spectral 
model above. To subtract the nuclear contribution, we used a local region at the same radii 
but with different azimuthal angles than the clouds. The results are listed in Table 2. Fig- 
ur eE] shows th e se [OII I] to soft X-ray ratios, comparin g the measurements w ith HRC results 



m 



Wang et al.l J2010bh and photoionization models in iBianchi et al.l ( l2006r . 



The [OIII] /X-ra y ratios at the two X-ray hot spot locations (corresponding to radio 
knots C2 and C5 in iMundell et J1 Il995|)_ were also nieasure d, and found to be uniformly 
low, ~3. This was previously noted in IWang et al.l (l2009bl ). and implies enhanced X-ray 
emission compared to other clouds under nuclear photoionization, which is likely associated 
with the jet-cloud interaction and consistent with their association with high NelX/OVII 
ratio regions. 
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2.4. Spatial Morphology of the Fluorescent Line Emission 

The 6.4 keV fluorescent iron Ka emission of neutral or mildly ionized cold material is 



readily visible in the ACIS spectra of the nucleus (IWang et al.ll2010bl ) and the extended region 
(Figure [7^). Another fluorescent feature, Sil Ka emission at 1.74 keV, is also clearly present 
in the extended emission (Figure Ub)- Both spectra of extended emission were extracted 
from a 2" < r < 30" annular region centered on the nucleus. Is this "extended" fluorescent 
line emission truly originating from outside the nuclear region, or is it the result of the PSF 
wing spreading out point-like nuclear features? 

To determi ne the spatial distributi on of the Fe Ka line emitting material, we followed 



the procedure in lSmith fc WilsonI ( 120011 ): a continuum image, taken to be the average of the 
counts in the 5.9-6.2 keV and 6.5-6.8 keV bands, was subtracted from the image extracted 
in the 6.2-6.5 keV band, resulting in a "pure" FeKa line image, shown in Figure |8^. The 
Fe Ka emission appears mostly circular in the r < 2", but shows faint extended emission 
towards the north and west. 



The Si Ka line is blended with the MgXII Ly/3 and SiXIII /, r emission ( lOgle et al. 



2000), thus the above continuum subtraction is not applied. We extracted emission in the 
1.7-1.8 keV range to obtain the Si Ka emission map (Figure (Hb)- It is mostly concentrated 
within a r ~ 2" circle, with a slight elongation along the optical bi-cone direction (northeast- 
southwest). 

However, at such energies the wings of the energy-dependent telescope PSF spread a 
small but non- negligible amount of the nuclear photons to large radii. Thus the presence of 
"extended" Fe Ka emission or Si Ka emission in the r > 1" region may in fact come from 
the PSF-scattered unresolved nuclear emission. 

To investigate this point, we further compared the observed radial profile of the Fe Ka 
emission band (6.2-6.5 keV) and the radial profi le in the same band for a point-like nucleus 
simulated with MARX0 (see Wang et al. 2010b), shown in Figure IHt- A similar comparison 



for the Si Ka (1.7-1.8 keV) is shown in Figure Eb. The extracted observed profiles shown in 
Figure [9] agree well with the PSF profiles, except for the inner r < 2" where the simulation 
overpredicted the observed emission due to heavy pile up in the PSF core. 

The simulated PSFs indicate that the r < 2" region contains 89% of the point source 
emission in the 6.2-6.5 keV band and 94% of the point source emission in the 1.7-1.8 keV 
band. Approximately 6% of the nuclear Fe Ka emission is expected to be spread within the 



''Available at 



http : //space . mit . edu/CXC/MARX/ 
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2" ^ ^ 6" extended regions, the extent of Fe Ka seen in Figure [7^. Taking into account 
the contribution from the PSF wings, the remaining extended Fe Ka emission is 5 ± 1% of 
the total Fe Ka emission. This should be considered an upper limit for the truly extended 
emission, because the nuclear Fe Ka emission is likely to be underestimated due to pile 
up. Similarly, we estimated that the extended Si Ka emission is < 4 ± 1% of the total Si 
Ka emission. As pointed out by the referee, the measured extended emission could easily 
disappear with a slightly broader simulated PSF given that the DitherBlur parameter in the 
MARX simulation is not known a priori. 

Therefore, we conclude that although there is faint extended emission seen in the Fe 
and Si fluorescent line images, at most 5% of the observed Fe Ka emission and 4% of the 
Si Ka emission can be truly spatially extended; most of these emissions are consistent with 
the PSF scattering of a strong nuclear component. 



3. Discussion 
3.1. Jet— Cloud Interaction in a Photoionized Outflow 



Both our previous work using HRC image of NGC 4151 nuclear region (jWang et al. 



2009bl ) and this work using ACIS find that most of the NLR clouds in the central 150 pc 
radius region of NGC 4151 have a relatively constant [OIII] to soft X-ray ratio (~10; Figure[6] 
and Table 2), no matter the distance of the clouds from the nucleus. This ratio indicates a 
uniform ionization parameter and a density decreasing as r~^. This may be consistent with 
a photoionized wind scenario (e.g., Bianchi et al. 2006) as the NLR clouds are outflowing 
(e.g.. Kaiser et al. 2000). The [OIII] /X-ray ratios at the locations of jet-cloud collision are 
lower (~3), leading to the suggestion of enhanced X-ray emission in these regions due to 
shock heating in addition to photoionization. 

This conclusion is supported by the results of multiwavelength imaging studies. The high 
spatial resolution emission line images of OVII, OVIII, and NelX emission show extended 
structures, in particular two "hot spots" in the NelX image, which are close to the rad i o knot 



features (Figures [2] and [3]). At these locations (1" from the nucleus), iMundell et al.l (|2003[ ) 
noted that the morphology of the highly collimated radio outflows appears more disturbed, 
suggesting interaction with some dense clouds here. There are also clearly optical cloud s 



with high velocity dispersion associated with the jet-cloud impa ct (Mundell et al. 2003[). 



The radio knots are near to the location of strong [Fell] emission ( IStorchi-Bergmann et al. 



20091 ). falling between the peaks of the [Fell] emission. 



Altogether, the morphology suggests a scenario where the radio ejecta runs into a 
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dense r cloud in the inhomogeneous ISM and results in locally enhanced shock heating 
(e.g., ICapetti et al.lll999l ). This jet-cloud interaction has been suggested to explain the 



enhanced [Fell] emission, which could be due to the iron unlocked from grains by the shocks 
( IStorchi-Bergmann et al.ll2009l ). The NelX emission, which traces the higher ionization gas, 
beco mes more prominent with the extra collisional ionization from the jet-cloud interac- 



tion. lArmentrout et al.l (120071 ) also noted that the poor fit to the OVII line profile in the 
grating spectra could be due to presence of a non-photoionization component. This possi- 
bility is further supported by our spectral fits, which suggest that the NelX/OVII enhance- 
ments have a thermal origin, requiring the presence of a collisionally ionized plasma (the 
kT = 0.58 ±0.05 keV com ponent) in addition to a photoionized nuclear outflow (lYang et al. 



20011 : IWang et al.ll2009bl l 



The emission measure of the APEC component allows us to estimate the electron 
density (~ uh) and the thermal pressure pth (~ 2nekT). The emitting volume (V^) of 
the hot gas, assuming that the depth along the line of sight is comparable to the other 
dimensions, is 5.6 x 10^*^ cm^ for the high NelX/OVII ratio hot spots. The filling factor rj 
is assumed to be 100% and both of rig and pth have a weak dependence on i] (p oc r]~^^'^). 
There is ~ 10% uncertainty in kT, which is propagated to the estimated pressure and energy. 
For the kT = 0.58 keV emission, we derived an average density rie = 0.06 cm~^, a hot gas 
pressure pth = 6.8 x 10"^" dyne cm~^, a thermal energy content of Eth = 4 x 10^^ ergs, and 
a cooling time of Tc = Eth/ Lapec ~ 10^ yr. The estimated internal pressure of the radio jet 
based on the synchrotron minimum energy assumption ranges between 10"'' and 10~^ dyne 
cm~^ (Pedlar et al. 1993), to be confined by the hot gas pressure of 10~^ dyne cm~^, which 
further supports our collisional ionization scenario. 



^Assuming T = 10^ K for the photoionized clouds and rig = 220 cm ^ (iPenston et al. 

1990l ). the thermal pressure from the photoionized gas is Pph = 3 x 10^^*^ dyne cm~^, which is 



comparable to the thermal pressure of the hot ISM, implying a possible pressure equilibrium 
between the collisionally ionized hot gas and the photoionized line-emitting cool clouds. 

We further estimate the age of the int eraction from the ap proximate crossing time of the 
0.5" region (dimension of the radio knot; iPedlar et al.lll993l ) with a characteristic velocity 
of Cs = 200 km s~^, the local sound speed. This is approximately the thermal velocity 
for the 0.58 keV gas (vt h = 240 km s~^), and a l so at the order of velocity dispersion of 
emission line gas seen in IStorchi-Bergmann et al.l (120101 ). Assuming the strong shock jump 
conditions, a crude estimate of the sh ock velocity Vgh can be obtained using T 1.5 x 10^ K 
(vsh/lOO km s-^f jRaga et al.ll2002l ). For the T ~ 6.7 x 10^ K X-ray-emitting gas, a Vsh of 
~700 km s~^ relative to the downstream material is required. This appears consistent with 
a supersonic but relatively slow radio jet, which is constrained to be sub-relativistic from 
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radio proper motion estimates (e.g., vjet < 12000 km s ^; lUlvestad et al.ll2005[ ). Thus we 
obtain a characteristic timescale of tcross ~ 10^ yr. 

If the jet-cloud interaction converts kinematic energy into heating of the hot gas, a lower 



limit can be placed on the kinematic luminosity of the jet, L^.e. ^ Eth/tc 



1 X lO'^^ erg 



s ^. It is interesting to compare this observed Lk.e. to the jet power directly inferred from the 



synch rotron emission. Using the Pjet^Pradio scaling relation (Equation 1) in lCavagnolo et al. 



( 20101) and a total radio luminosity at 1.4 GHz {uL,^ 



4 X 10^'' erg s ^; 



Mundell et al. 



19991 ) ■ we find Pjet ~ 10^^ erg s ^, suggesting that > 0.1% of the jet power is deposited in 



the ISM. However, recalling that the P jet^Pmdin relation was d erived from a sample of radio 
loud galaxies with X-ray cavities (see ICavagnolo et al.l l2010l for review on the Pjet^Pradio 
scaling relation) while here NGC 4151 is of low radio power, w e also estimate directly the jet 
energy flux using the pressure in the knots ( iPedlar et al.lll993l ) based on the minimal energy 
assumption, giving a similar Pjet = 5 x 10^^ erg s~^. 

We further note that, the crossing time tcross is comparable to the cooling time Tc, 
implying that the hot gas is prominent only locally (close to the radio knots). Indeed, along 
the path of the radio jet, the clouds that are spatially close to the jet impact spots (e.g., clouds 
#3,4 bracketing knot C2; #9,10 for knot C5) show [OIII]/X-ray ratios indistinguishable from 
others (Figure [6]). If these clouds are physically related to the ISM interacting with the jet, 
instead of being projected close to the knots in our line-of-sight, their X-ray emission shows 
little evidence for excess shock heating in addition to nuclear photoionization. This localized 



heati ng is perhaps related to the highly collimated radio ejecta in NGC 4151 flMundell et al. 



2003), in contrast to the more expanded lobe-like radio outflow in Mrk 3 (jCapetti et al. 



19991 ) ■ where the NLR gas appears broadly impacted. We plan to investigate in future work 
whether the [OHI] /X -ray ratio, such as <4 measured here and seen in Mr k 3 and NGC 1386 
(jBianchi et al.l l2006l ) . together with other line ratios (e.g., [FeII]/[PII]; lOliva et al.l 1200 ll ) 
could be useful diagnostics for strong jet-ISM interaction in Seyfert galaxies. 



In terms of energetics, we summarize our findings in NGC 4151: 
estimated from the radio power is 

7 x 10^3 



(a) The jet power 
10^^ erg s~^, a few percent of the current AGN en- 



ergy output (L 



bol 



erg s 



-1. 



Kaspi et al.l |2005| ) , which is close to the canonical 
5% of the radiant energy adopted in the oretical models for efficient AGN feedback (e.g. 



Pi Matteo et al.lbood : faopkins et al.lbooeh : (b) > 0.1% of the jet power is observed to have 
been deposited into the host ISM in the nuclear region through the interaction between the 
radio jet and the dense medium; (c) The shock heating due to the jet-cloud interaction ap- 
pears localized to the impact spots, and most clouds are consistent with being photoionized 
by the nucleus. 
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3.2. Extent and Origin of the Fe Ka Line 



Ogle et al.l ( 120001 ) claimed that the narrow iron line emission was spatially resolved in 



the Chandra HETG observation and 65% ± 9% of the Fel Ka emission comes from the 
ENLR at distance up to 6"extent (~400 pc across ). This is intrig uing, since the putative 
pc-scale torus or more compact broad-line region (ILiu et al.ll2010l ar id references t h erein) 
are expected to be the primary location of narrow Pel Ka emission. ISchurch et al.l (120031) 



Ogle et al. (l2000h 



cautioned that due to the sensitivity of the short grating observation in 
the 1-3 arcsec region off the peak of cross-dispersion profile contains few photons. Therefore 
the detection of spatially extended iron line emission cannot be significant. A high signal-to- 
noise XMM-Newton spectrum of the neutral Fe Ka suggests that all lin e flux originates in a 
nearly Compton-thick tor us that is not resolvable at current resolution (ISchurch &: Warwick 
2OO2I : ISchurch et al.ll2003l ). 



Our results show that the extended Fe Ka emission, if present, is only at ~ 5% level 
of the observed Fe Ka emission, when the contribution from the PSF wings is taken into 
account. This is consistent with our constraint on the extent of the Si fluorescent emission 



and supports the conclus ions of ISchurch et al.l (120031 ) . while it is in st r ong d isagreement 
with the 65% reported in lOgle et al.l J2OO0I ). We note that in lOde et al.l ( l200oh the PSF is 
represented as a narrow Gaussian with FWHM = 0.9", when the Chandra mission was new. 
The MARX model is now much more advanced so that a more reliable estimate becomes 
possible. At 6.4 keV, the encircled energy fraction of a point source is only 63.7% (see 
Chandra Proposers' Observatory Guid^ Chapter 4 Table 4.2) in such a PSF core. Hence 
a significant fraction of the nuclear emission could have been measured in Ogle et al. as 
spatially extended Fe Ka emission. 



4. Conclusions 

In this paper we present spectral analysis and emission line images from deep Chandra 
observation of NGC 4151, aiming to resolve and characterise the X-ray emission in the inner 
~130 pc-radius nuclear region. The findings are summarized as follows: 

1. We have obtained high spatial resolution X-ray narrow-band images of OVll, OVllI, 
and NelX line emission, which are blended at the ACIS spectral resolution. The 
images show extended structures that are spatially correlated with the radio outflow 



^http : //cxc ■ harvard . edu/proposer/PDG/| 
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and optical [OIII] emission. 

2. We find strong evidence for jet-ISM interaction, including morpliological correspon- 
dences witli regions of X-ray enliancement, peaks of NIR [Fell] emission, and optical 
clouds. This is further strengthened by the presence of a kT = 0.58 keV collisionally 
ionized component in the spectral fitting of the hot spots, and by the excess of X- 
ray emission in addition to nuclear photoionization as indicated by a low [OIII] /X-ray 
ratio. We find a possible pressure equilibrium between the collisionally ionized hot 
gas and the photoionized cool clouds. The estimated velocity of the shocks from the 
jet-cloud impact is ~700 km s~^. 

3. We estimate that the jet power in NGC 4151 is close to a few percent of the current 
AGN energy output {Lboi ~ 7 x 10^^ erg s^^). The derived thermal energy in the 
hot gas suggests that > 0.1% of the jet power is deposited into the host ISM in the 
nuclear region through the interaction between the radio-jet and the dense medium. 
The [OIII] /X-ray ratios of NLR clouds bracketing the radio knots show little deviation 
from other photoionized clouds, indicating a localized impact on the ISM by the highly 
coUimated jet. 

4. We investigate the spatial extent of the fluorescent features, including the Fe Ka emis- 
sion and the Si Ka emission. Our results show that both are dominated by point-like 
emission, consistent with an origin of unresolved inner structure such as a torus. The 
extended Fe Ka emission is < 5% of t he observed Fe Ka emission, which is in strong 



disagreement with the 65% reported in lOgle et al.l ( 120001 ). 
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Fig. 1. — The tri-color composite image of the central 7" x 7" (450 pc-across) circumnuclear 
region of NGC 4151, where the soft (0.3-0.7 keV), medium (0.7-1 keV), and hard band (1-2 
keV) images have been smoothed with a FWHM = 0.3" Gaussian kernel and shown in red, 
green, and blue, respectively. The contours are of the medium band and highlight a linear 
feature spatially coincident with the radio jet. The pixel scale is 0.0625" per pixel, 1/8 the 
native ACIS pixel. 
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Fig. 2. — The emission- line structure of the hot ISM in the central ~250 pc of NGC 4151 (a) 
OVII; (b) OVIII+OVII RRC; (c) NelX. The nucleus position is indicated with a cross. The 
contours outline the radio outflow in the 1.4 GHz MERLIN map (Mundell et al. 1995). The 
ACIS images have been rebinned to 0.0625" per pixel, and smoothed with a FWHM = 0.3" 
Gaussian kernel. 
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Fig. 3. — Same emission-line maps of the central ~250 pc of NGC 4151 as in Figure [21 but 
with contours showing the H ST /FOC [OIII]A5007 emission line clouds (Winge et al. 1997). 
(a) OVII; (b) OVIII+OVII RRC; (c) NelX. The ACIS images have been rebinned to 0.0625" 
per pixel, and smoothed with a FWHM = 0.3" Gaussian kernel. 
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Fig. 4. — (a) A close view of the NelX emission line image overlaid with the contours of 
VLB A+VLA radio jet (green lines Mundell et al. 2003) and NIR [Fell] (black dashed 
lines; IStorchi-Bergmann et al.l (120091 )). Note how the NelX line emission peaks appears to 
align with the general direction of the linear radio outflow, (b) Ratio image between the 
NelX band and the OVII band images. Contours are the same as in (a). The arrows (red) 
indicate the knots of high NelX/OVII ratios, and the circles (magenta) mark the spectra 
extraction regions. The box (magenta) marks the extraction region for the low NelX/OVII 
ratio region. The crosses (yellow) mark the locations of the HST clouds with high velocity 
dispersion (Kaiser et al. 2000; see also Figure 7 in Mundell et al. 2003). 
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Fig. 5. — Spectral fit of tlie X-ray emission in liot spot (black) and surrounding regions (red). 
Li kely identificatioii s of tlie blended emission lines are labeled, based on the HETG line list 
in lOgle et al.l ( 2000 ). See text and Table 1 for model details. 
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Fig. 6. — The [OIII] to soft X-ray ratio as a function of the cloud's distance to the nucleus. 
The triangles are the measurements in Wang et al. (2009), and the squares are from this 
work. The filled (red) circles represent the two high NelX/OVII ratio regions (Figure H]). 
The blue, red, green, and magenta lines indicate the [OIII]/X-ray ratios for NGC 7212, NGC 
3393, Mrk 3, and NGC 1386 (Bianchi et al. 2006), respectively. The dotted, dashed, solid, 
and dot-dashed lines are the Cloudy model predicted values from Bianchi et al. (2006) for 
different radial density profiles: rif. cx where /3 = —1.6, —1.8, —2, and —2.2, respectively. 



5.5 6 6.5 1.6 1.7 1.8 1.9 

Energy (keV) Energy (keV) 

Fig. 7. — (a) The spectrum of the extended emission showing presence of Fe Ka line, (b) 
Same spectrum but showing the Si Ka hne part. The underlying sohd hne is the NGC 4151 
nuclear spectrum (modeled with a power law plus a gaussian), scaled to 6% and 3% of the 
nuclear flux, respectively. 
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Fig. 8. — (a) The contour map of continuum-subtracted Fe Ka line emission of the central 
15" X 15" region in NGC 4151. (b) The contour map of Sil Ka line emission of the same 
region, extracted between 1.7 keV and 1.8 keV. The cross marks position of the nucleus. The 
contour levels are of logarithmic scale, between 3 x 10~^° to 4 x 10~^ photons cm~^ s"^. 
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Fig. 9. — (a) Comparison between the radial profile of simulated PSF in the 6.2-6.5 keV 
band (red dashed line) and the observed narrow-band Fe Ka line emission (black solid line), 
(b) The same comparison but for the 1.7-1.8 keV Si Ka line emission band. 



Table 1. Spectral Fits to the NelX/OVII Enhancements and Surrounding Region 



Region 


Counts 


logC/i 


\ogNH,i 


Norml 


logC/2 


log Nh, 2 


Norm2 


kT 


Norm3 


Nuclear PSF=' 






[0.3-2 kcV] 














[keV] 






[erg s"! cm— ^ arcscc— ^] 


High-ratio 


4585 


-0.01 ±0.03 


20.5 


8.0 X 10-1^ 


1.9 ±0.2 


22.5 


1.5 X 10-1* 


0.58 ±0.05 


1.8 X 10-5 


1455 ± 35 


8.9 X 10-1"' 


Low-ratio 


1317 


-0.01 ±0.03 


20.5 


8.2 X lO-i** 


1.9 ±0.2 


22.5 


6.3 X 10-19 


0.58 ±0.05 


4.5 X 10-^ 


530 ± 23 


6.8 X 10-14 



^The absorption column Njj is fixed at 2 X 10 cm" -2 jMurphv et al.lll996l l. the Galactic column towards NGC 4151; / d.o.f.= 146/142 for the best fit presented here. 
''Expected counts from simulation of PSF scattered nuclear emission in the 0.3-2 keV band. 
'^Soft X-ray surface brightness of the extended emission after nuclear emission removed. 



Table 2. Measured X-ray and [OIII] Fluxes 



Cloud 


Distance 


Distance 


[OIII] flux 


Net Counts 


0.5-2 kcV Flux'' 


[OIII] /soft X 


[OIII] /soft X 


label 


to Nuc. (") 


to Nuc. (pc) 


(10~^^erg s~l cm~^) 


(0.3-2 kcV) 


(lO-l^'erg s-l cm-2) 


(ACIS) 


(HRC) 


1 


1.76 


114 


1.1 


369±24 


1.5 


7.4 


12 


2 


1.44 


93.6 


1.9 


688±29 


2.7 


6.9 


12 


3 


1.01 


65.6 


2.2 


581±38 


2.3 


9.4 


12 


4 


0.78 


50.7 


1.3 


470±39 


1.9 


6.9 


7 


5 


1.0 


65.0 


1.5 


168±29 


0.7 


22.2 


75 


6 


0.6 


39. 


1.1 


99±44 


0.4 


27.8 


110 


7 


0.4 


26. 


3.4 


1320±96 


5.3 


6.4 


5 


8 


0.4 


26. 


3.3 


560itl00 


2.2 


14.7 


13 


9 


0.85 


55.2 


2.4 


694±47 


2.8 


8.6 


3 


10 


0.89 


57.8 


1.6 


266±35 


1.1 


15.0 


8 


11 


2.15 


139.7 


0.6 


264±17 


1.1 


5.7 


15 


C2'' 


0.9 


58.5 


0.8 


1180±44 


4.7 


1.7 


2 




0.93 


60.4 


1.4 


1052±43 


4.2 


3.3 


3 



''An average conversion factor 7.2 X 10 erg s ^ cm (counts 8 ^) ^ is derived from spectral models in Table 1 to estimate flux from net count rate. 

''C2 and 05 are the radio knots in Mundell et al. (1995), corresponding to the locations of high NelX/OVII ratio hot spots to the east and west of the nucleus, 
respectively. 



